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ABSTRACT 


We  construct  in  this  paper  a general  purpose  algorithm  for  sol- 
ving polynomial  0-1  programming  problems.  The  algorithm  is  applied 
directly  to  the  polynomial  problem  in  its  original  form.  Further, 
no  additional  variables  are  introduced  in  the  solution  process. 

The  algorithm  was  tested  on  randomly  generated  modest  size  prob- 
lems and  the  preliminary  computational  results  obtained  are  very 
encouraging. 


1. 


INTRODUCTION 


We  consider  in  this  paper  the  general  (linear  or  nonlinear)  0-1 
programming  (GP)  problem  of  the  form 


T 

Max  c x 


(GP) 


s.t. 


f±(x)  < b±  (i  = 1,  ...  , m) 

x.  € (0,1)  (j  = •••  , n) 

J 


where  the  f.'s  are  polynomials  in  the  x.  variables. 

1 J 

It  was  shown  that  every  (GP)  problem  is  equivalent  to  a special 
structured  linear  0-1  problem  - the  Generalized  Covering  (GC)  problem 
of  the  form 


T- 

Min  c x 


(GC) 


s • <«*  • 


Ax  > b 
1 


x.  = x e (0,1} 
J J 


(j=l,  ...  , n) 


where  A = (a^ ^ ) with  e {0,l,-l}  and  b^  = 1 + ]^MLn{a.^,0} 

J 

(i  = 1,  ...  , m).  See,  for  example,  Granot  and  Hammer  [5,6]  for  the 
general  (GP)  problem  and  Balas  and  Jeroslow  [1]  for  the  linear  (GP) 
problem. 

Granot  and  Hammer  [5]  attempted  to  employ  this  equivalence  relation 

in  an  algorithm  for  solving  (GP)  problems.  Unfortunately,  this  attempt 

failed  because  of  the  large  number  of  constraints  in  the  equivalent 

n n 

(GC)  problem.  For  example,  a constraint  of  the  form  x.  < [-k]  is 

j=l  J 

equivalent  to  f n 'i  generalized  covering  constraints. 

I 
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Another  approach  for  solving  nonlinear  (GP)  problems  suggests  to 
first  linearize  the  nonlinear  constraints  and  then  solve  the  equivalent 
linear  problem,  see  e.g.  [2,3,8].  However,  the  most  severe  shortcoming 
of  this  linearization  approach  stems  from  the  radical  increase  in  the 
dimension  of  the  problem.  For  example,  a (GP)  problem  with  n 0-1 
variables,  m constraints  and  k distinct  products  of  variables  is 
transferred,  by  Glover  and  Woolsey  [3J>  to  an  equivalent  linear  problem 
with  n + k variables  and  m + n + k constraints.  Clearly,  this 
approach  will  produce  unmanageable  problems  for  large  k. 

A branch  and  bound  algorithm  for  solving  linear  and  nonlinear  0-1 
problems  was  devised  by  Hammer  in  [7]* 

In  [4]  an  efficient  general  purpose  algorithm  (not  a branch  and 
bound  type)  was  developed  for  solving  the  monotone  (GP)  problem,  i.e., 
when  the  coefficients  of  the  f^'s  are  all  nonnegative.  In  this  paper 
we  extend  the  results  of  [4]  and  construct  an  algorithm  for  solving 
general  (GP)  problems.  The  algorithm  solves  a sequence  of  relatively 
small  (GC)  problems,  where  each  (GC)  problem  is  a relaxation  of  the  original 
(GP)  problem  as  well  as  a tighter  relaxation  than  its  predecessor  in  the 
sequence. 

Though  the  algorithm  is  applicable  to  both  linear  and  nonlinear 
(GP)  problems,  its  promise  lies  mainly  for  the  nonlinear  case.  One  of 
its  main  advantages  is  that  no  additional  variables  are  introduced  in 
the  solution  procedure. 

Preliminary  computational  results  for  randomly  generated  (GP)  problems 


are  reported  in  Section  4. 


2.  PRELIMINARIES 


rr 


* 


Consider  first  the  linear  (GP)  problem  of  the  form: 


(1) 


T 

Max  c x 


s.t. 


(2) 

O) 


Ax  < b 
x e (0,1} 


, ^nxl  . mxl  . mxn  . ...  _ „ 

where  ceR  , beR  , AeR  are  given  with  c > 0. 


Every  constraint  of  (2)  can  be  replaced  by  an  equivalent  set  of 
generalized  covering  constraints  in  the  following  manner.  Consider 
a single  inequality  of  the  form: 


(M 


J = 1 


where  |a'|  > |a^|  > •••  > |a^|.  For  every  j for  which  a^  < 0 we 


can  substitute  x.  = 1 
J 


x to  obtain  an  equivalent  constraint  of  the 

J 


form 


(5) 


n a. 


a . 


Z aixi,]  = Z lal!xiJ  < b'  - Z min(0,a  ) = b 
j=l  J J j=l  J J j=l  J 


where 


a. 
„ J 


xj  lf  aj  - 0 


if  a3  < 0 * 


A subset  VCN  = {l,2,  , n}  is  said  to  be  a cover  of  (5)  if 


(6) 


jeV 


aj>b 


It  is  said  to  be  a prime  cover  if  no  proper  subset  of  it  is  a cover. 
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n a. 

If  we  further  denote  by  S = [x:  £ a.x.^  < b}  and  by  ft  the  set 

j=l  J d S 

of  all  prime  covers  of  (5),  then  it  is  easy  to  see  (e.g.  [1,5]) 

that 


a,  _a. 

(7)  x e S iff  7T  x.J  = 0 V V.  c iff  £ x.J  > 1 V V.  e 0 


JeV. 


_a . a . 

where  x .^  = 1 - x . . 

Applying  the  above  procedure  to  every  constraint  in  (2)  will  result 
with  an  equivalent  generalized  covering  problem. 

This  result  can  be  further  generalized  (see  [5])  to  (GP)  problems 
of  the  form 


T 

Max  c x 


(8) 


s.t. 


f^(x)  < K (i  = 1,  m) 


Xj  e (0,1) 


(j  “'ll  ...  i n) 


where  f^x^,  •••  > xn)  are  polynomials  of  the  form: 


(9) 


[i<v  - > V *ik x 


k=l  jeN. 


ik 


with  any  subset  of  N. 

In  order  to  obtain  the  set  of  generalized  covering  constraints 
which  is  equivalent  to  a polynomial  inequality  of  the  form 


(10) 


= I aV  w , xj  1 b’  > 


k=l  jeN. 


one  can  proceed  as  follows:  For  every  k for  which  a£  < 0 substitute 
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£ x.  7T  x.  for  TT  x.  in  (10)  to  obtain  an  equivalent  constraint 

jeN£  J J 

l < j 


of  the  form: 


£ a TT  x j a'  7r  x + £ |a'|  £ x.  TT 


k=l  JeN, 


k jeN*  J ieN’ 

■«  < J 


< b*  - £a*  = b 


where  Z (Z  ) denotes  summation  over  positive  (negative)  values  of  a^ 
+ 
a. 

and  x.J  equals  to  either  x.  or  x.  in  accordance  with  (11).  Now, 

J J *3 

ai 

by  substituting  y = V x.  in  (ll)  we  obtain  an  equivalent  linear 

jc»k  ■> 

inequality  in  the  variables  y^  of  the  form: 


Z \ yk  1 b • 


Using  the  method  described  earlier  one  can  now  obtain  the  set  fl 

b 

of  prime  covers  to  (12).  Then  (12)  is  equivalent  to 

a 

V y = 0 V V e fi  iff  V IT  x J = 0 V V.  € . 

keVi  K 5 keVi  jeNk  3 

Thus,  an  inequality  of  the  form  (lo)  is  equivalent  to  the  following  set 
of  generalized  covering  inequalities 


_a 

x 3 > 1 V V.  € i7 
j€Vi  J * 1 S 


where  V.  = N N.  . 
1 keVj. 


Employing  the  above  procedure  to  every  constraint  of  (GP)  results 
with  an  equivalent  generalized  covering  problem. 

3.  AN  ALGORITHM  FOR  SOLVING  (GP)  PROBLEMS 

We  construct  in  this  section  a general  purpose  algorithm  for  solving 
(GP)  problems.  We  first  generate  a generalized  covering  problem  which 
is  a relaxation  of  (GP),  referred  to  as  the  Generalized  Covering  Relaxation 
(GCR)  problem.  If  the  optimal  solution  to  (GCR)  is  feasible  to  (GP), 
it  is  also  optimal.  Otherwise  (GCP)  is  augmented  by  additional  generalized 
covering  constraints  that  eliminate  this  optimal  solution.  Each  augmented 
constraint  will  be  shown  to  dominate  the  corresponding  generalized  covering 
constraint  produced  by  the  method  of  Section  2.  The  (GCR)  problem  is 
then  resolved.  The  variables  in  each  (GCR)  problem  are  the  original 
variables  of  (GP)  and  its  constraints  are  a small  subset  of  the  constraints 
of  the  equivalent  (GC)  problem. 

A distinctive  feature  of  this  algorithm  is  that  it  does  not  generate 
the  entire  equivalent  (GC)  problem.  Rather,  the  size  of  each  (GCR) 
problem  solved  is  relatively  small.  Additional  constraints  are  generated 
only  to  eliminate  the  optimal  solution  to  (GCR)  when  not  feasible  to 
(GP).  Moreover,  no  additional  variables  are  introduced  in  the  solution 
procedure.  The  variables  in  each  (GC)  problem  are  the  original  variables 
of  the  (GP)  problem. 

Consider  again  the  general  (GP)  problem: 


T 

Max  c x 

s.t. 
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and  assume,  without  loss  of  generality,  that  c.  > 0.  Let  (GCR)  be  any 

J 

relaxation  of  (l4)  with  x its  optimal  solution.  Further  let  I denote 
the  set  of  indices  corresponding  to  constraints  in  (l4)  which  are  violated 
by  x.  For  each  i e I denote  by: 

st(x)  = (k : a > 0 and  ir  x.  = 1} 

• «T  J 


S7(x)  = {k:  a < 0 and  TT  x = 1} 

1 J 

S?(x)  = {k:  a < 0 and  V x.  = 0}  , 

1 1 

i.e.,  st(x)(s7(x))  is  the  set  of  indices  corresponding  to  terms  in 

the  i 1 constraint  with  positive  (negative)  coefficients  that  do  not 

vanish  at  x,  and  S?(x)  corresponds  to  terms  in  the  ith  constraint 

with  negative  coefficients  that  do  vanish  at  x.  Let  S^(x)  denote 

the  ordered  set  S^(x)  = (s(l),s(2),  ...  , s(£)}  where  s(i)  e sT(x)  U 3?(x) 

and  s(k)  < s(j)  if  either  |ais^|  > lais(j)l>  or  if  lais(k)l  = 


lais(j)l  and  k<j. 


Consider  now  the  constraint 


u aisCk)  1 
k=l  ls(k;  jcN. 


is(k) 


jeSi(x) 


which  was  derived  from  the  i constraint  in  (l4),  i e I#  Since  (15) 
is  violated  by  x,  it  can  serve  for  generating  generalized  covering 
constraints  which  when  augmented  to  (GCR)  will  eliminate  x.  The 
constraint  we  will  generate,  using  the  following  procedure,  dominates 
that  constraint  produced  from  an  original  violated  constraint  using  Granot 
and  Hammer's  method  [5]  (as  described  in  Section  2). 
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Observe  that  since  Procedure  I is  applied  only  to  constraints  i, 
i e I,  the  procedure  will  terminate  after  applying  Step  2 at  most 
|s^(x) | times,  where  |s^(X)|  is  the  cardinality  of  S^(x).  Moreover 
(l6)  dominates  the  corresponding  constraint  obtained  using  Granot  and 
Hammer's  method  [5].  This  since  a term  in  (15)  with  a negative  coefficient 
contributes  at  most  a single  variable  to  (l6). 


I 
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Example  1: 

Consider  the  following  constraint: 

Jjx-jXg  - Jx^Xj^x^  + 2x^Xg  < 4 

and  the  solution  x = (1,1,1, 1,0,0) . Since  x violates  the  given 
constraint,  we  apply  Procedure  I to  the  following: 

(IT)  5x^2  - Jx^x^  < 4 

which  will  result  with  the  generalized  covering  constraint: 

(18)  + x2  - x^  > 0 . 

Applying  the  method  described  in  Section  2 for  generating  a general- 
ized cover,  we  first  transfer  (17)  to  the  following  equivalent  constraint 
with  positive  coefficients 

5x1x2  + 5x3  + 5x3x4  + 5x3x4x5  1 7 * 

Now  since  x^  = x^  = 0 we  will  derive  a generalized  cover  from 

5x^2  + 3x^X5  < 7 

which  results  with 

Xi  + x2+X3+X4-x5>0. 

This  constraint  is  obviously  weaker  than  (l8). 

Lemma  1: 

Augmenting  (OCR)  with  (16)  eliminates  x and  results  with  a new 


(GCR)  problem  which  is  a relaxation  of  (GP) 
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Proof: 


"fch 

Since  x violates  the  i constraint  in  (1*0,  we  have 
P, 


(19) 


b < £ aik  n 

jeNik 


x . = > a . . + a.  . 

3 + ^ i3  r *3 

jeSi(x)  jeSi(x) 


But 


i 


(20)  5,  a«  * S a^<*>  j./  * b 

jeS^(x)  k"1  JeNis(k) 


jeS7(x) 


a.  . = b. 
13  i 


l ' ' — 1’ 

and  thus  x violates  (15).  Any  cover  derived  from  (15)  will  result 
with  a constraint  that  when  augmented  to  (GCR)  will  eliminate  x. 
Further,  since  (l6)  is  valid  for  (15),  it  is  also  valid  for  (l4)  and 
thus  no  feasible  solutions  to  (l4)  are  eliminated. 


An  Algorithm  for  Solving  (GP)  Problems: 

Step  0:  Start  with  (GCR)  being  the  trivial  relaxation,  i.e., 

Max{c  x:  x c {0,l}),  whose  optimal  solution  is  x = (1,  ...  , l). 

Step  1:  Let  I be  the  set  of  indices  corresponding  to  constraints  of 
(14)  which  are  violated  by  x.  Apply  procedure  I to  generate 
k generalized  covering  constraints  for  each  i e I and  add 
them  to  (GCR). 

Step  2:  Solve  the  augmented  (GCR)  problem  to  obtain  its  optimal  solution 
x.  If  x*  = 1 - x^  (j  = 1,  ...  , n)  is  feasible  to  (GP), 
terminate  with  x an  optimal  solution  to  (GP).  Otherwise 


go  to  Step  1. 
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Theorem: 

The  algorithm  described  above  converges  to  an  optimal  solution  of 
(GP)  in  finitely  many  iterations. 

Proof: 

From  Lemma  1 we  conclude  that  (GCR)  is  a relaxation  of  (GP),  which 
implies  optimality.  In  step  1 of  the  algorithm  we  eliminate  an  optimal 
solution  to  (GCR)  which  is  not  feasible  to  (GP).  Since  the  number  of 
binary  vectors  is  finite,  the  convergence  follows . 


Remark  1: 

Observe  that  (15)  is  violated  by  x for  every  i e I and  thus  every 
one  of  the  k X |l|  (where  J I | is  the  cardinality  of  i)  constraints 
added  to  (GCR)  is  sufficient  to  eliminate  x.  Generating  more  than  one 
generalized  cover  will  hopefully  reduce  the  number  of  (GCR)  problems 
needed  to  be  solved . 


Remark  2; 

Note  that  no  additional  variables  are  introduced  by  the  algorithm 
in  the  solution  process.  Further,  no  special  devices  are  needed  to  handle 
nonlinear  constraints. 


Example  2: 


Solve 
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Max  Jx^  + 4x2  + 6xj  + x^  + x^  + 5xg 
s.t. 

7x1  - 2x2x4  + ^x3  - 2x1+  - x5  + 6xg  <4 

4xxXj  - 8x2Xj  + Jx^gX^  - 2x^  _ Jx^Xg  + 7x^  < -1 

3Xix2  + 4xj  - x^x  < 5 

5x^x^  + 4x2xg  - 4x.jX^  + 2Xj  - x^x^xg  < 4 

x^,  x2,  Xy  X^,  x,.,  Xg  € (0,l)  . 

We  start  with  x = (1,1,1, 1,1,1) . 

Iteration  I: 

For  x = (1,1, 1,1, 1,1)  the  left-hand  sides  equal  (12,1,6,6), 
hence  every  constraint  is  initially  violated,  i.e.,  I = {l,2,J,4). 

We  start  by  setting  k = 1 and  thus  generate  one  generalized  covering 
constraint  from  every  violated  constraint.  Since  sT(x)  = 4>,  i e I, 
the  generalized  covers  are  derived  from  the  following  constraints: 

7x1  + 4x^  + 6xg  < 9 

4x1Xj  + 3x1x2Xj  + 7xg  < 12 

5x1X2  + ^X3  - 6 

5x1Xj  + 4xgX^  + 2Xj  < 4 . 

The  (GCR)  problem  obtained  is 
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Min  3x1  + 4xg  + 6x^  + x^  + x ^ + Sxg 

s.t. 

*1  + x6  > 1 

*1  + *3  + *6  -1 

X1  + x2  + ^3  > 1 

*1  +X2  + X^  + Xg  > 1 

x.^  e {0,1}  (i  = 1,  ...  ,6) 

whose  optimal  solution  is  x^  = 1,  x.  = 0,  (i  = 2,  ...  ,6). 

Iteration  2: 

For  x = (0,1,1, 1,1,1)  the  left-hand  sides  equal  to  (5, -6, 3, 5). 
Since  the  only  violated  constraints  are  the  first  and  the  last  one, 

I = {1,4}  and  we  derive  additional  generalized  covers  from  the  following 
constraints 

Ux3  + 6x6  - 9 

4x2x6  " 4x1x4  + 2x3  - 5 * 

When  augmenting  the  newly  derived  generalized  covers  to  the  previous 
(GCR)  problem  we  obtain,  after  some  simplifications, 

Min  3x^  + 4xg  + 6xj  + x^  + x^  + 5xg 

s.t. 

X1  + x6  > 1 

x3+x6  > 1 

xl  + x2  + x 3 > 1 

x2  + x3  + x6  - xu  > 0 

Xj^  € {0,1}  (i  = 1;  •••  , 


6) 
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whose  optimal  solution  is  x.^  = x(  = 1,  x2  = x^  - x^  = x,-  = 0.  Since 
x = (0,1,1, 1,1,0)  is  feasible  to  the  original  problem,  it  is  also  an 
optimal  solution. 

Observe  that  if  we  attempt  to  solve  the  above  example  by  first 
transferring  it  to  a linear  0-1  problem,  via  Glover  and  Woolsey' s most 
efficient  transformation  (3],  the  resulting  equivalent  problem  will 
have  11  additional  variables  and  17  additional  constraints. 

4.  SIMPLIFICATION  AND  COMPUTATIONAL  RESULTS 

It  is  well  known  that  simplification  and  elimination  rules  can  be 
applied  to  the  rows  and  columns  of  the  coefficient  matrix  of  a covering 
problem  before  solving  the  problem.  Similar  simplification  and  elimi- 
nation rules  can  be  devised  for  the  generalized  covering  problem.  These 
simplification  rules,  to  be  described  below,  are  employed  to  any  newly 
generated  generalized  cover  in  order  to  reduce  the  size  of  the  coefficient 
matrix  A of  the  (GCR)  problem. 

Let  a^(a^)  denote  the  i^*1  row  column)  of  the  matrix  A.  Then 

the  following  simplifications  are  applied  at  each  iteration,  see  also  [5]. 


Simplification  1 

"fch 

If  a row  a^.  contains  a unique  1 (unique  -l)  in  the  j column, 

then  set  x.  = 0 (x.  = 1)  in  the  original  (GP)  problem. 

J J 


Simplification  2 

If  there  exist  rows  a^,  a^  for  which  aik  ^ 0 implies  aik  = a..^, 

then  row  a.  can  be  eliminated. 

J 


Simplification  3 

If  there  exist  rows  a^,  a^ 
aii  = 1 = "aj i>  then  substitute 


for  which  a..  = a..  V k / 1 and 
ik  ,ik 


a . = 0 and  eliminate  row  a.  from  A. 
jl  1 
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Simplification  4 

If  there  exist  rows  a.,  a. 

and  aik  ^ 0 ^P11*53  aik  = ajk’ 

and  eliminate  row  a^. 


and  index  p such  that 
k / p,  then  substitute 


iP 


Simplification  5 

If  there  exist  rows  a.,  a.,  and  indices  r,  p for  which  a = 

i j iK 

a..  = 0,  k / r,  p and  a = -a  . = -a.  = a.  ^ 0,  then  substitute 

jk  lr  jr  ip  JP 

xr  = x^  in  the  original  (GP)  problem. 


Simplification  6 

If  there  exist  rows  a.,  a.  and  indices  r, 

1 J 

a.,  =0,  k i r,  p and  a.  = -a . = a.  = -a. 

jk  * r > y ir  jr  ip  jp 

x = 1 - x in  the  original  (GP)  problem, 
r p 


p for  which  = 

^ 0,  then  substitute 


Simplification  7 

If  there  exist  a row  a.  for  which  a..=0Vk^i,  r and 


*U  ■ air  ' 1 

Ki  = air  = 

-D, 

then  in  any  row  a^  for  which 

‘ -Jr  * 1 

substitute 

ajr  = 

0 (ajjg  = 0), 

and  in  any  row  a . 

J 

for  which  a ^ ^ 

= -a,  = -1 

jr 

substitute  a^  = 0 

<v  - o)- 

Simplification  8 


If  there  exist  a row  a.  for  which  a.,  = 0,  k / l,  r and 

1 IK 

a . . = -a.  = 1,  then  in  any  row  a.  for  which  a..  = a.  =1  (a., 

i !■  ir  J J*  Jr  j-e 

a . = -l)  substitute  a . 

jr  jr 


0 (a^.O). 


Simplification  9 

If  there  exist  columns  ad,  for  which  a^  > a*^  and  < c^, 

then  column  a^  can  be  eliminated  (i.e.,  x^  = 0 in  the  current  (GCR) 
problem) . 


Simplification  10 


If  there  exist  a column  a and  an  index  set  J 


for  which 


then  column  a can  be  eliminated 


Our  algorithm  was  programmed  in  Fortran  IV,  implemented  on  an 
IBM  370/l68  and  tested  on  randomly  generated  problems. 


In  the  preliminary  computational  results  to  be  reported,  the  vector 


uniform  on  [0,10].  Every  constraint  is  of  the  form 


a. .T.  < b.  with 


the  number  of  terms  is  chosen  uniformly  between  3 and  K 


the  number  of  variables  in  each  term  is  chosen  uniformly 


is  chosen  uniformly  between  -20  and  20 


is  chosen  uniformly  between  1 and  the  number  of  variables 


summation  over  negative  (positive)  coefficients.  Thus 
[ Li»Ui ] denotes  the  range  of  the  left-hand  side  of  the 
ith  constraint. 


a constant 


the  number  of  constraints 


the  number  of  variables 


By  changing  a,  we  tested  the  influence  of  large  and  small  values  of 


b on  the  computation  time 


IT  - 


In  the  following  four  tables  we  summarize  the  preliminary  compu- 
tational results.  Each  cell  is  averaged  from  4 runs  and  contains 
three  entries,  A,  B,  C,  which  designate 

A - Number  of  iterations  (i.e.,  (GCR)  problems  solved) 

B - Number  of  covers  in  the  last  (GCR)  problem  solved  (after  sim- 
plifications) 

C - True  execution  time  in  seconds,  i.e.,  total  CPU  time  minus 
loading,  input/output  and  overhead  time. 


a = .6 

a = .75 

a = .9 

11,  28,  1.30 

* 

* 

li,  32,  2.1*3 

3.6,  9-3,  .07 

* 

12,  34,  4.68 

6.6,  15.5,  .11 

1,  3.6,  .03 

Table  1;  M=20,  K = 7,  b^  = Gfth  + (l-oOl^  (i  = 1,...,20) 
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— 

; m = 15 



M = 20 

M = 25 

0 

r<V 

II 

X 

1 5,  13,  .16 

7,  20,  .32 

9.5,  32.5,  1.14 

14,  54,  4.4  : 
1 

Table  3:  N = 4o,  K = 7i  b*  is  uniform  on 

^■12  Ui  + 12  ^i*  + J2  ^ i ^ = • • • Jl) 


I 

j M = 15,  N = 30 

M = 20,  N = 40 

*1 

M = 25,  N = 50 

6.2,  17,  .25 

10,  26.5,  2.14 

— ■■  - ■ ■ ■ . t 

17,  36,  6.38 

4-  - — * 

Table  4 : K = 7#  b+  is  uniform  on 

- ~ + 22  ^2  + 22"  ^ = • • • > M) 

The  following  conclusions  can  be  drawn  from  the  preliminary  compu- 
tational results. 

a)  Tables  1 and  2 reveal  that  the  tighter  the  constraints  are  the 
more  difficult  it  is  to  solve  a (OP)  problem.  We  recall  that  the  covering 
relaxation  algorithm  for  solving  monotone  0-1  problems  [4]  has  precisely 
the  same  property. 

We  remark  that  when  we  attempted  to  further  restrict  the  r.h.s. 
vector  b (i.e.,  decreasing  a)  almost  all  of  the  randomly  generated 
problems  were  found  to  be  infeasible  by  our  algorithm.  Further,  increasing 
the  number  of  constraints  M beyond  the  ratio  n/m  = 2 also  results 
with  almost  all  randomly  generated  problems  being  infeasible. 

The  r.h.s.  b^  in  Tables  3*  4 was  chosen  so  as  to  roughly  vary  between 
2/3  and  3A  of  the  range  R^. 
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b)  By  comparing  Tables  1 and  2 we  conclude  that  the  difficulty 
of  solving  a (GP)  problem  increases  with  an  increase  in  the  number  of 
terms  in  each  constraint  (i.e.,  an  increase  in  k).  This  result  again 
was  expected  since  by  increasing  K we  increase  the  number  of  equivalent 
generalized  covering  constraints. 

c)  Tables  1 and  2 exhibit  also  the  significant  effect  of  an  increase 
in  the  number  of  variables  on  the  execution  time.  We  remark  though  that 
the  increase  in  both  the  number  of  iterations  and  the  number  of  generalized 
covers  is  more  moderate. 

d)  Table  3 reveals  that  an  increase  in  the  number  of  constraints 
increases  significantly  the  execution  time.  Again,  the  increase  in  the 
number  of  iterations  and  the  size  of  the  (GCR)  problems  is  more  moderate. 

e)  Table  4 exhibits  the  effect  of  increasing  both  M and  N,  in 
a fixed  ratio,  on  the  solution  characteristics.  Again,  the  increase 
in  execution  time  is  very  significant  while  the  increase  in  the  number 
of  iterations  and  the  number  of  generalized  covers  is  more  moderate. 

It  appears  that  the  genralized  covering  relaxation  approach  for 
solving  modest  size  (GP)  problems  is  quite  promising.  For  example, 

(GP)  problems  with  50  variables,  20  constraints  and  K = 10  art  solved, 
on  the  average,  in  12.81  seconds.  What  is  particularly  encouraging  is 
the  relatively  modest  number  of  iterations  required  to  solve  a (GP) 
problem,  for  example , ( GP ) problems  with  N = 50  , M = 20  , K = 10  are 
solved  on  the  average  in  20  iterations.  Clearly,  by  improving  the  ef- 
ficiency of  our  generalized  covering  algorithm  we  can  further  reduce 
the  execution  time. 
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Some  improvements  of  the  algorithm  are  being  currently  considered. 
In  particular,  devising  and  employing  approximate  algorithms  for  both 
the  ( GP ) and  the  (GCR)  problems  in  order  to  accelerate  the  implicit 
enumeration  algorithm  used  for  solving  the  (GCR)  problem. 
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